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We conduct a systematic investigation of the nuclear collective dynamics that emerges in systems 
with random two-body interactions. We explore the development of the mean field and study its 
geometry. We investigate multipole collectivities in the many-body spectra and their dependence 
on the underlying two-body interaction Hamiltonian. The quadrupole-quadrupole interaction com- 
ponent appears to be dynamically dominating in two-body random ensembles. This quadrupole 
coherence leads to rotational spectral features and thus suggests the formation of the deformed 
mean-field of a specific geometry. 



PACS numbers: 21.60.Cs, 21.60.Ev, 24.60.Lz 
I. INTRODUCTION 



Emergent phenomena is one of the most profound top- 
ics in modern science addressing the ways that collectiv- 
ities and complex patterns appear from multiplicity of 
components and simple interactions. Ensembles of ran- 
dom Hamiltonians allow one to explore the emergent phe- 
nomena in a statistical way, and thus to establish generic 
relations and rules. To study the many-body physics of 
interest we adopt a shell model approach with a two-body 
interaction Hamiltonian. The sets of the two-body inter- 
action strengths are selected at random resulting in the 
two-body random ensemble (TBRE). Symmetries, such 
as rotational, isospin, and parity, entangled with complex 
many-body dynamics result in surprising regularities dis- 
covered recently in the low- lying spectrum. Patterns ex- 
hibited by the random ensembles are remarkably similar 
to those observed in real nuclei. The high probability for 
the ground state spin to be zero is the most astounding 
feature of the TBRE discovered in Ref. [I] . Signs of al- 
most every collective feature seen in nuclei, namely, pair- 
ing superconductivity, deformation, and vibration, have 
been observed in random ensembles f2^'6j . While the sys- 
tematics of the ground state quantum numbers is almost 
not sensitive to the short-range pairing matrix elements, 
the probability to find a coherent paired structure in the 
wave- functions of low-lying states is enhanced [2]. The 
presence of rotational features in the spectra is another 
unexpected result seen in the TBRE [51 [7j. 

The goal of this work is to study the emergence of col- 
lective mean-field dynamics in ensembles with random 
interactions. The discussion is organized as follows: In 
Sec. |IT]we briefly deflne the TBRE, introduce signatures 
of collective motion, and discuss ways to detect them. In 
Sec. |III| we present our study of collectivities in singlej 
level models. More complex models are explored in Sees. 
llVl andlVl We summarize our results in Sec lVII with a dis- 
cussion of the quadrupole-quadrupole Hamiltonian which 
appears to be responsible for most of the observed phe- 
nomena. 



II. COLLECTIVE OBSERVABLES AND 
MODELS 

In the spirit of the traditional shell model approach, 
we deflne a model conflguration as (ji, j2 . . . where 
N nucleons occupy a set of single particle levels labeled 
by their angular momentum j. In this work we assume 
that the single particle energies are degenerate. We ex- 
amined other models for which this was not the case and 
the results are similar. The Hamiltonians in the TBRE 
are defined with a set of two-body matrix elements which 
are selected at random. The distribution of the matrix 
elements is Gaussian so that, within a given symmetry 
class, the ensemble of Hamiltonian matrices for two parti- 
cles coincides with Gaussian Orthogonal Ensemble. The 
presence of rotational symmetry and, where relevant, of 
parity and isospin symmetries is assumed. The typical 
number of random realizations was between 10^ and 10^ 
for all ensembles presented in this work. 

In the TBRE the number of realizations where the 
ground state spin Jgg = is disproportionally large. 
Aiming at collective phenomena we select realizations 
with Jgs = 0. With the exception of the ground state, la- 
beled as Ogs, we denote the low-lying states by the value 
of their spin with an identifying subscript. The subscript 
is given in bold font if it refers to the absolute position 
of a given state in the spectrum. Throughout the paper 
we give probabilities of finding realizations with certain 
features, these probabilities are always quoted in percent 
relative to the size of the ensemble; however, all proba- 
bility distribution plots are normalized to unit area. 

In order to identify and to analyze manifestations of 
collective phenomena in the spectra we use a set of ob- 
servables. The goal is to choose a finite number of spec- 
tral observables that are likely to convey the most in- 
formation about possible collective structures in a scale- 
independent way and with minimal model dependence. 
These quantities and the logic behind their selection are 
discussed in what follows. 

The geometry of the nuclear mean field is described by 
the multipole density operators Mx^ with multipolarity 
A and magnetic component /i. The structure of the mul- 
tipole operators depends on the valence space, for each 
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model it is addressed separately. The reduced transition 
rate from an initial state \JM) to a final state \J' M') 
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is one of the observables. Here \JM) denotes a many- 
body state with angular momentum J and magnetic pro- 
jection M. The total transition strength from a state J 
is given by the sum rule 

5a(J) = ^ B{EX, J^J') = {JM\Y,M\^Mx,\JM), 

.7' M 
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which provides a convenient normalization to assess the 
fractional collectivity of the transition 
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The shape of a state is described by its multipole mo- 
ments specified by the expectation value 



Qx{J) = (JJ|Xao|JJ>. 
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For a non-spherical system this moment describes the 
shape of a deformed nucleus measured in the lab frame. 
The intrinsic shape is characterized by the body-fixed 
(intrinsic) multipole moments Q\. A rotational spectrum 
(band) emerges for every fixed intrinsic shape. In a rigid 
rotor these intrinsic moments are the same for all states 
in the band and they determine the lab- frame observables 
in Eqs. ([ij and (|4]). For the ground state band of inter- 
est, the intrinsic moments determine the total transition 
strength Sx{Qgs) = Ql- 

In the axially symmetric case the quantum number K, 
a projection of the angular momentum onto the body- 
fixed symmetry axis, is conserved. Then for each rota- 
tional if -band the relations between the observables in 
the lab frame and in the intrinsic frame are expressed via 

- - xa.jj^^ 
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Clebsch-Gordan coefficients Qx{J) = Q^^'io jj^xo jk 



This limit of an 



and B{E\,J J') = 

axially symmetric rotor provides a convenient normal- 
ization to examine the multipole moments. In this work 
instead of Qx{J) we quote a normalized intrinsic moment 
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which is computed as if the state is a member of the 
K = rotational ground state band. 

In this paper we only briefiy touch the subject of col- 
lectivities other than quadrupole, see Sec. |III C[ thus 
for convenience the subscript A is omitted for A = 2. 
The relation between the lab- frame moment of the 2i 
state and its intrinsic moment is Q(2i) = — 2/7Q(2i). 
For the axially symmetric rotor the quadrupole transi- 
tion sum rule for the Ogs is saturated by a single tran- 
sition b{E2,0gg — > 2i) = 1. The quadrupole moment is 
q(2i) — 1 for prolate or q{2i) — —1 for oblate shapes. 



We normalize the total transition strength Sx{J) to 
its maximum possible value for a given valence space. 
Taking the A = 2 case as an example, we define the 
quadrupole-quadrupole (QQ) Hamiltonian as 



(6) 



The eigenstate energy of the QQ Hamiltonian coincides 
with the total transition strength ^ for that state: 
^QQiJ) — —S{J). Thus, the absolute value of the 
ground state energy of the QQ Hamiltonian \EQQ{Ogs)\ 
is the maximum possible value of the total transition 
strength S{J) for a given model space and for a given 
structure of the quadrupole operator. We therefore de- 
fine a relative transition strength as 
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To summarize, in our study we use the dimensionless 
variables defined in Eqs. ([s]), (|5|, and Q. To shorten 
notations we define b = b{E2,0gs — t- 2i), q = q(2i), and 
s = s(Ogs). For collective models of pairing, rotations, 
and vibrations 5 ~ 1. We refer to a realization with 
b > 0.7 as collective and with b < 0.3 as non- collective. 
The quadrupole moment q allows one to separate differ- 
ent collective modes: q ss ±1 for rotations and g w for 
vibrations and for paired states. In what follows we al- 
lude to collective realizations with q > 0.7 as prolate and 
those with q < —0.7 as oblate. For rotations the rela- 
tive transition strength s is proportional to the square of 
the intrinsic moment, and thus it is associated with the 
Hill- Wheeler deformation parameter Within Elliot's 
SU(3) model [8^ the relative transition strength s can be 
thought to represent the expectation value of the Casimir 
operator which identifies the irreducible representation. 
In cases where s sa 1 the ground state band structure is 
close to that of the QQ Hamiltonian. 

The collective structure is further analyzed using the 
following 4i state. The types of collective modes can be 
classified by the ratio of the excitation energies measured 
relative to the energy of the Og^ state 
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This ratio is close to for pairing, 2 for vibration, and 
10/3 for rotation. The ratio of deexcitation rates 
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B{E2Ai^'2i) 
B{E2,2i^0gs) 



(9) 



is another measure. It is nearly for pairing, 2 for vibra- 
tional mode, and 10/7 for rotational motion. Typically, 
for models with the QQ Hamiltonian R42 and ,642 are 
close to the rotational values, see summary in Tab. |T] 
A comprehensive review of different collective models, 
their analytic predictions, and comparisons with rota- 
tional spectra observed in real nuclei can be found in the 
textbooks 
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III. THE SINGLE j LEVEL MODEL 
A. Quadrupole collectivity 

We begin our presentation with single j level mod- 
els. Starting from the original paper [l^ the single j 
level with identical nucleons has been at the center of 
numerous investigations; a good summary may be found 
in the following reviews |51 HJ |5J [TT]. With many is- 
sues understood and with still unanswered questions, the 
single j model remains an important exploratory bench- 
mark. The model, while simple, has a number of partic- 
ularly attractive features which can be of both advantage 
and disadvantage [12 ■. the Hamiltonian is defined with a 
small number of parameters; apart from an overall nor- 
malization constant, the multipole operators are uniquely 
defined; a special role is played by the quasispin SU(2) 
group; and the particle-hole symmetry is exact. 

In Fig. [T] the system with 6 nucleons in a single 
j 19/2 level is examined, we refer to this system as 
(19/2)^. Here we select 10.4% of random realizations 
where the Ogs state is followed by the 2i state. The distri- 
bution of the fractional collectivity b = b{E2,0gs — > 2i) 
in Fig. [l|^a) points to highly collective nature of the 
quadrupole transition Ogg 2i. Most realizations with 
Ogs and 2i are collective (6 > 0.7), their fraction is 7.8% 
of the total number of samples. These realizations are 
shaded in red. This collectivity is not a statistical coinci- 
dence. The system (19/2)^ has 1242 spin-states, among 
them there are 10 states with J = and 23 states with 
J = 2. Thus, statistically the chance for the O^s, 2i spin 
sequence to occur among all other possible outcomes is 
only 0.015%. The large fractional collectivity for the 
transition between these two states is even more unlikely, 
given that the transition strength is shared among 23 
J = 2 states, the chances for b{E2,0gs — ?► 2i) > 0.7 are 
of the order of 1 in 10^. 

There are two peaks in the distribution of the 
quadrupole moment in Fig. [ijb), they reflect prolate 
and oblate deformations. For most of the collective re- 
alizations, which are shaded in Fig. [T| the magnitude of 
the quadrupole moment is consistent with the value for 
the axially deformed rigid rotor {\q\ w 1). The ground 
state is most likely to be oblate, but in about one out of 
four collective cases a prolate mean field emerges. 

The collective realizations are further analyzed in 
Fig. |2] where the distribution of the relative transition 
strength s is shown. In Fig. [2] the oblate (q < —0.7) 
and prolate {q > 0.7) cases are shaded with different 
patterns. The relative transition rate s for the oblate 
samples is close to the maximum possible s = 1. Thus, 
for these realizations the ground state band structure is 
similar to that of the QQ Hamiltonian. The data on the 
QQ Hamiltonian for our models is summarized in Sec |VI[ 
For prolate systems the distribution of the relative tran- 
sition strength peaks around s = 0.37. 

In Fig. [3] we focus on distributions of the deexcitation 
ratio ,642 and the energy ratio R24 defined in Eqs. [9] 
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Figure 1: (Color online) (19/2)*^. (a) The distribution of the 
fractional collectivity b. (b) The distribution the intrinsic 
quadrupole moment q. Only realizations with the Ogs, 2i spin 
sequence are included in both panels. There are 10.4% of such 
realizations. The 7.8% of collective realizations (6 > 0.7) are 
shaded. 

and [Sj We use the same shading for prolate and oblate 
realizations as in Fig. |2] but slightly modify our selection 
of samples. We chose the collective realizations that have 
states 2i and 4i, with 2i being not higher than the forth 
excited state and 4i state being above it. 

The collective oblate realizations comprise a peak in 
the distribution of B42 in Fig- [sf^a) around the rota- 
tional limit of B42 — 1.4. For prolate realizations the 
distribution peaks near ,642 = 0.8 and has an extended 
shoulder. It is likely that the rotational structure is frag- 
mented in instances with a weak prolate deformation. 
Here the lower value of s seen in Fig. |2]is used to suggest 
a weak deformation. Therefore, the 4i state is not purely 
rotational. 

The distribution of the ratio of the excitation energies 
-R42 in Fig. |3|b) seems to contradict the rotational limit. 
For most of the collective realizations the values of the 
ratio fall between the pairing limit of 1 and the vibra- 
tional limit of 2, while in the rotational limit the ratio 
of 3.3 is expected. This discrepancy has been reconciled 
in Ref. p] with the observation that the rotational or- 
dering emerges for the ensemble-averaged excitation en- 
ergies. The same conclusion is expected from the geo- 
metrical chaoticity arguments [13 . Excitation energies 
are sensitive to non-collective features, this leads to large 
fiuctuations of i?42. The experimental observations of 
realistic nuclei also show that when the quadrupole tran- 
sition rates follow the rotational systematics, the exci- 
tation energy spectrum can deviate from rotational; on 
occasions, the spectrum is closer to the vibrational limit 
[S]. The coexistence of both prolate and oblate config- 
urations in this (19/2)^ system could be another reason 
for the distortion in the energy spectrum. Within Elliot's 
SU(3) model analogous mixing of group representations 
was investigated in Ref. [14 . 



B. Triaxiality 

The triaxiality is marked by the presence of low- 
lying levels 22, 3i, 42, 5i. The excitation energies 
are subject to equalities E{2i) + E{22) = £^(3i) and 
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Figure 2: (Color online) (19/2)''. The distribution of the 
relative transition strength s for the collective realizations 
(shaded area in Fig. [T|. The quadrupole moments shown 
in the inset are separated into prolate {q > 0.7) and oblate 
{q < —0.7) shapes. The resulting distributions are shaded 
with a pattern and a uniform color, respectively. The frac- 
tion of oblate cases is 5.2% and the fraction of prolate cases 
is 1.3% relative to the total number of random realizations. 
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Figure 3: (Color onhne) (19/2)**. (a) The distribution of the 
deexcitation ratio £42 defined in Eq. ||9|. (b) The distribution 
of the excitation energy ratio _R42 defined in Eq. ([8|. The 
distributions are comprised of 13.6% of realizations that have 
the Oga, 2i, 4i sequence with b > 0.7, the 2i state is not 
higher than the fourth excited state, and -E(4i) > _E(2i). The 
prolate cases and oblate cases, that appear in the ensemble 
with probabilities 3.3% and 7.1% respectively, are shaded with 
the same patterns as in Fig. |2] The values of 542 and i?42 
for the QQ Hamiltonian listed in Tab. [I] are marked with the 
vertical grid lines. 



This angle is determined by the three reciprocal moments 
of inertia Ai, i = 1,2,3 in the rotor Hamiltonian 



tan2r 
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The ratio of the excitation energies of the 2\ and 22 states 
is another parameter of the rotor Hamiltonian. It is con- 
venient to express this ratio i?2i22 iii terms of the angle 
7£)i? defined using the Davydov-Filippov model of irro- 
tational flow jTS] as 



sin2(37DF) 



2(l+i?2i2. 
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In our example that follows, the triaxiality is small and 
Idf ~ 0.5i?2i22- Thus, the rotor Hamiltonian, given 
by the three moments of inertia, can be equivalently de- 
scribed by an overall energy scale, the iC- mixing angle F, 
and the angle 7df- 

The quadrupole shape is parametrized by the Hill- 
Wheeler parameters /3 and 7 which deflne the quadrupole 
operator M.2^i- The relation between the parameters of 
the rotor Hamiltonian and the intrinsic shape is model- 
dependent. The irrotational-flow moments of inertia 
discussed in Ref. [15] result in 7£)i? — 7 and F — 
{arccot [8001(87)] ~ 7} /2; the latter implies F ^ 7 for 
small triaxiality. A rather different result follows from 
the rigid-body moments of inertia. 

We determine F, 7Di?, and 7 independently from the 
spectroscopic observables. The parameter ^-^p is ob- 
tained from the energy spectrum, Eq. (10). Following 
Ref. [TH] one can view the sum rules 



2i) + 6(^2,0g, ^22) = 1 



and 
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6(£;2,2i^22)+g2(2i) = l 



for the J = 2 two-state model as the Pythagorean the- 
orem for amplitudes. The angles in the corresponding 
right-angled triangles are 7 — F and 7 + 2F, therefore 



tan2(7-F) = 



5(^2,0^22) 
B{E2,0^ 2i)' 
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4:E{2i) + E{22) = E{5i). It is remarkable that these re- 
lations appear to be well satisfied by the spectrum of the 
QQ Hamiltonian in the (19/2)^ configuration, for which 
i?2i3i + R2^3^ = 1.005 and 4i?2i5i + R2,5, = 1-026. Here 
Rjj' = E{J)/E{J') denotes the ratio of excitation en- 
ergies. The rigid rotor Hamiltonian, defined by three 
moments of inertia, is responsible for these correlations 
in the spectrum. 

In this work we examine two low-lying 2i and 22 states. 
These are the only states with spin 2 in the triaxial rotor 
model, they are mixed configurations of if = and K = 
2. We use angle F to express the level of the if- mixing. 



tan2(7 + 2F) = 



2B{E2,2i 22) 

WW) ■ 
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These equations allow one to determine the triaxiality 7 
and the iC- mixing angle F. 

All three angles 7, F, and ^df are small in our models 
with the QQ Hamiltonian, see discussion in Sec. |VI[ 
Correspondingly, in the TBRE the effects of triaxiality 
are weak but detectable. 

For our studies of triaxiality presented in Fig. |4]we use 
the (19/2)^ model. We recall that in the triaxial rotor 
model there is a second 22 state with Q(2i) = — Q(22). 
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Thus, we select collective realizations with Ogs and 2i, 
and in addition to that we require that in the entire 
spectrum there is a 22 state for which the equality 
2(22) = — Q(2i) holds within 20% of accuracy. In collec- 
tive realizations of rotational type the magnitude of the 
(5(22) is large as compared to the quadrupole moments 
of other many-body states. This simplifies the identifica- 
tion of the 22 state. We find that practically for all collec- 
tive realizations this second 22 state exists. Indeed, from 
the total number of random realizations a large fraction, 
18.3%, satisfy all of the mentioned triaxiality conditions. 
In Figs. |4][a) , |4|^b) , and |4|^c) we show the distributions 
of the triaxiality angle 7, if-mixing angle F, and 7df, 
respectively. We use the same shading as in Fig. [3] to 
separate prolate and oblate shapes. 

In the (19/2)^ model one often finds collective realiza- 
tions with oblate intrinsic deformation and s ~ 1, these 
realizations are triaxial with 7 « 9°, Fig. ga). This 
result, as well as F sa 0° in Fig. |4|b), is consistent with 
that of the QQ Hamiltonian. The less frequent prolate 
cases are nearly axially symmetric. 

In the TBRE the angle '^df, Fig. |4j^c), appears on 
average to be higher than the corresponding angle in the 
QQ Hamiltonian. The peak in the ^df distribution is 
also higher than the peak in the 7 distribution, compare 
Figs. |4ja) and[4|c). (We remind that ■foF = 7 in the ir- 
rotational flow model.) Nevertheless, no conclusions can 
be made from these two discrepancies. We believe that 
the excitation energies could be influenced signiflcantly 
by non-collective features. The situation may be similar 
to the one in Fig. |6jb), where 4i state is lower than ex- 
pected for the rotor. Similarly, if the 22 state is lowered 
the resulting ^df is larger. In both cases the lowering is 
relative to the excitation energy of the 2i state. 



C. Higher multipole moments 
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It is known that in the TBRE the probability to 
flnd a Ogs state followed by either one of the states 2i, 
4i, 61, or 81 is disproportionally large as compared to 
what is statistically expected. For the (19/2)^ model 
the corresponding probabilities are 10.4%, 17.3%, 11.9% 
and 1.8%. In an attempt to understand this, we re- 
peat the previous study but target the collective real- 
izations of multipolarity A = 4, 6, and 8. For realiza- 
tions with the Ogs state and with the first excited state 
of spin A in Fig. [5] we consider the fractional collectiv- 
ity &(A) = b{E\,Ons ~> Ai) and the multipole moment 
q\ = q\{Xi). Fig. p] shows evidences for intrinsic shapes 
with deformations of higher multipolarities. In partic- 
ular, for A = 4 and 8 there is a sizable number of col- 
lective realizations where 5(A) > 0.7. These realizations 
are shaded (in red). The corresponding distributions of 
the multipole moments in Figs, [sj^d) andjsjf) have peaks 
which are centered at non-zero values of qx. The A = 6 
shape collectivity is nearly absent in the (19/2)^ system: 
the realizations are mostly non-collective, 6(6) < 0.7 



Figure 4: (Color online) (19/2)'^. (a) The distribution of the 
triaxiality angle 7 . (b) The distribution of the if-mixing 
angle F. (c) The distribution of the triaxiality angle yoF 
from the Davydov-Filippov model. The angles are obtained 
from Eqs. (Ill, (12 1 and (10 1. We select realizations with 



two states of spin 2 in the spectrum and require b > 0.7 and 
q(2i) ~ —5(22); 18.3% of realizations satisfy this set of re- 
strictions. The realizations with prolate and oblate shapes 
are shaded with the same patterns as in Figs. [2] and [3] Ver- 
tical grid lines indicate the triaxiality parameters calculated 
from the QQ Hamiltonian, which are: 7 = 9.79, F = 0.73, 
and 7I5F = 7.52. 



(shaded in blue), and the corresponding moment has a 
peak centered near zero. 

Investigations of other single j systems show presence 
of multipole collectivities with A = 2, 4, 6, and 8. Gen- 
erally, the collectivities corresponding to the intrinsic 
quadrupole shape are the most pronounced ones, however 
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Figure 5: (Color online) (19/2)*^. The distributions of the 
fractional collectivity 6(A) are shown in panels (a), (b), and 
(c). The distributions of the intrinsic multipole moments 
qx are shown in panels (d), (e), and (f). The plots are or- 
ganized in three rows corresponding to multipolarities with 
A = 4, 6 and 8. Here we include realizations where, in ad- 
dition to the Ogs state, the first excited state is either 4i, or 
6i , or 8i . The shaded areas correspond to collective and non- 
collective modes with fe(A) > 0.7 and 6(A) < 0.3 respectively. 
We use the same patterns as in Fig. [l] 

there are signatures of realizations with shapes of higher 
multipole deformation. The existence of such geometric 
structures may be related to the symmetries discussed in 
Ref. [U 



D. Multipole structure of the Hamiltonian 
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Figure 6: (Color online) (19/2)^. The same figure as Fig. 
[l] but for the 8-particle system, (a) The distribution of the 
fractional collectivity 6. (b) The distribution of the intrinsic 
quadrupole moment q. The histogram is comprised of 7.5% 
of random spectra with Ogs and 2i states. Shaded areas cor- 
respond to 4.6% of collective realizations (6 > 0.7) and 1.9% 
of non-collective realizations (6 < 0.7). This figure is analo- 
gous to Figs. [l]and[5] and the same shading is used in these 
figures. 




Figure 7: (Color online) (19/2)**. The same figure as Fig. [2] 
but for the 8-particle system. The distribution of the relative 
transition strength rule s for the collective realizations. This 
figure is analogous to Fig. [2] and the same shading is used as 
in Figs. [2] and [3] however only oblate shapes {q < —0.7) are 
seen. 



In this subsection we discuss the multipole structure of 
the two-body Hamiltonian in the single j level model. For 
this purpose we use a larger system of 8 nucleons in the 
same j = 19/2 model space, i.e. the (19/2)^ model. The 
distributions of the fractional collectivity, the quadrupole 
moment, and the relative transition strength shown in 
Figs. [6] and [7] are similar to the distributions observed in 
the (19/2)6 model in Figs. [l]and[2] The main difference 
between the models is that, in contrast to Fig. [ijb), 
only oblate ground state configurations are present in 
Fig. ^b). 

The collectivities observed in the single j studies are 
deeply rooted in the underlying geometric structure of 
the Hamiltonian. To focus on this relation we express 
the two-body Hamiltonian in the particle-hole channel in 
terms of the multipole operators 



^ = E^^E^L>'k- (13) 

K X 

The interaction parameters V)c in the particle-hole chan- 
nel are determined from those in the particle-particle 
channel Vz, via Pandya transformation 

Vk-^(2L + l)x£l4. (14) 

L 

The transformation coefficients 

K _ \ 3 j ^ \ 
^ 1 J J i J ' 
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are given by the six-j recoupling coefficients. 

On a single j level only even values of the particle-pair 
angular momenta L are allowed by the Fermi statistics. 
Thus, there are j + 1/2 interaction parameters Vl in Eq. 
( 14 1 . In studies of the TBRE a set of these parame- 
ters can be viewed £lS £l random vector in the j + 1/2 
dimensional space. There is no such a restriction on the 
particle-hole momentum IC. Thus, the inverse transfor- 
mation 



yL = ^(2/C + l)x£Vk 



(15) 



may produce some unphysical Vl with odd values of L. 
Such Pauli-forbidden terms in the Hamiltonian do not 
generate any dynamics. Therefore the 2j + 1 parameters 
Vjc contain passive components which can be removed 
making Vjc linearly dependent [12) . 

The interaction terms that correspond to the multi- 
poles with momentum JC = and K. = 1 are constants of 
motion [17 . The /C = term in Eq. (13 1, describes the 



nucleon-nucleon interaction that is the same for all an- 
gular momentum channels, Vl = x'l — const, as follows 
from Eq. (15). The resulting monopole Hamiltonian is 



proportional to the number of particle-pairs in a system. 
This Hamiltonian has no dynamical effect. Thus, there is 
no change in results if one constrains the TBRE by pro- 
jecting out the monopole JC ~ component as follows 



Vr ^ Vr 



.o J:L'i^L' + l)xl,VL' 
EL'i^L' + l){x'i,f 



(16) 



This effectively reduces the number of independent pa- 
rameters Vl- 

In a single j model space the /C = 1 multipoles are pro- 
portional to the angular momentum operators A4i^ ^ 
Jj^. Therefore the K. = 1 interaction leads to a rotational 
E{J) ^ J{J + 1) spectrum with Vic=i determining the 
moment of inertia. In the particle-particle channel, the 

operator is obtained with Vl = xi ~ const -I- L(L -1-1). 
Consistently, it was argued in Refs. j2][IH] that those in- 
teractions that lead to the positive moment of inertia are 
likely to result in the Jgg = 0. The exact operator 
component in the interaction can be removed by orthog- 
onalization to Xl following the procedure in Eq. (16). 

The changes in dynamics are no longer trivial when 
the quadrupole JC = 2 component in the interaction is 
modified. The role of different multipoles in the TBRE 
is studied in Fig. [8] and [9] whe re we remove different JC 
components from the interaction Hamiltonian in Eq. ( 13 ) 



using the Graham-Schmidt projection procedure. In the 
particle-particle channel the projection of pairing inter- 
action Vl = Sl,o has been extensively discussed in Ref. 
[2^ . The removal of pairing does not lead to any signifi- 
cant qualitative change, we thus forgo this topic in what 
follows. 

The probability to observe a certain ground state spin 
in the (19/2)* system is shown in Fig. [s] Three cases 
of random ensembles are reviewed: (a) the traditional 
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Figure 8: (Color online) (19/2)*. Probabilities to observe a 
certain ground state spin Jgs for three random ensembles: (a) 
the TBRE, (b) the TBRE without a term (the JC = 1 term 
in Eq. ([l3| is removed), and (c) the TBRE without both, 
1 and JC 



2 terms in Eq. ( 13 1 are 



and QQ terms (the JC 
removed) . 



TBRE where all j + 1/2 interaction parameters Vl are 
random Gaussian variables, (b) the case where JC = I 
term is removed, and (c) the ensemble where JC ~ 1 
and JC = 2 multipole components are removed from the 
Hamiltonian. While the wave functions in ensembles (a) 
and (b) are identical, the ground state spin distributions 
are different. The role of the moment-of-inertial-like 
term has been discussed before in Refs. |5] [131 [TS]; it 
appears to be fully responsible for the cases with max- 
imum possible spin. As seen in Fig. |8] the states with 
the maximum spin almost never occur as ground states 
in ensembles (b) and (c) where the interaction term 
(/C = 1) is removed. 

The ensembles (b) and (c) shown in Fig. [S] appears 
to have similar ground state spin distributions but the 
behavior of the fractional collectivity is different. In Fig. 
|9] for all three ensembles we show the distribution of 
the fractional collectivity of the transition between the 
Ogs and 2i states. It is evident that the quadrupole 
collectivity disappears once the quadrupole component 
in the interaction is removed. Thus, we conclude that 
the quadrupole-quadrupole component in the interaction 
generates the corresponding deformation and is respon- 
sible for the rotational behavior observed. 



IV. MODELS BEYOND SINGLE j 

In this section we expand the scope of our models 
and consider systems with two single-particle levels. The 
richer geometry allows one to study the effects of particle- 
hole conjugation, different structures of the multipole 
operators, and the role of parity. The distributions of 
the fractional collectivity and of the quadrupole mo- 
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Figure 10: (Color online) (13/2+, 13/2+)^ (a) The distri- 
bution of the fractional collectivity b. (b) The distribution 
of the intrinsic quadrupole moment q. The 4.1% of samples 
have the OgsiSi sequence. Shaded areas correspond to 1.2% 
of collective realization and to 1.8% of non-collective realiza- 
tions. This figure is analogous to Figs. [T] [5j and |6j and the 
same shading is used in these figures. 



Figure 9: (Color online) (19/2)*: The distribution of the frac- 
tional collectivity b for the same three random ensembles as 
in Fig. |8] Namely: (a) the traditional TBRE, (b) the TBRE 
without a term, and (c) the TBRE without both, and 
QQ terms. We select realizations with the Ogs state followed 
by the first excited state 2i. The fraction of such cases for 
ensembles (a), (b), and (c) is 7.6%, 8.2%, and 4.7% respec- 
tively. 




ment are shown in Figs. 10 for the (13/2+, 13/2"'")^ sys- 
tem. Here the model space is comprised of two levels 
with ji — j2 — 13/2. Both single-particle levels have 
the same positive parity, so that the effective spherical 
Hartree-Fock mean-field Hamiltonian can contain terms 
of a mixed structure such as aJ^Oj^. These terms are 
scalars for ji — j2. There is some arbitrariness in the 
choice of the single-particle matrix elements of the mul- 
tipole operator A^2p which depend on the radial overlap 
of the operator r^. We choose the radial overlap to be 



diagonal {ji\r^\ji) = {j2\r^\j2) 



and (ji 

,2 



= 0; other 

possibility with (jilr^lji) = (j2|r^|j2) = (jik^li2) has 
been explored and led to no substantial difference. 

A structurally different (13/2+, 13/2")^ model is ex- 
amined in Fig. [TT] where two levels of equal spin and 
different parity are considered. In this case the matrix 
elements of the Hamiltonian are restricted by parity. The 
same structure of the quadrupole operator is used. The 
model space of this kind has been explored in Ref. |T3] 
because it is the simplest model space that allows for 
quadrupole and octupole modes. The prevalence of the 
positive parity ground states is remarkable in this model. 
The ground state is most likely to have spin-parity 0+, 
4+, or 24+ with 35%, 19%, and 14% probabiHty, respec- 
tively. In contrast, 0~, the most probable negative parity 
ground state, happens only in 3% of realizations. For the 
(13/2+, 13/2~)^ model the number of many-body states 
is the same for both parities, 8,212 each. 

For both (13/2+,13/2+)6 and (13/2+, 13/2-)6 mod- 
els, the results related to the quadrupole collectivity are 
almost identical, see Figs. [TO] and [TT] Moreover, these 
results are similar to those for the single j level models. 



Figure 11: (Color online) (13/2+, 13/2")^ (a) The distribu- 
tion of the fractional collectivity 6. (b) The distribution of the 
intrinsic quadrupole moment q. The 6.9% of samples have the 
Ogs state and the first state 2i, both of positive parity. The 
2.4% of collective and 2.6% of non-collective realizations are 
shaded with patterns. This figure is analogous to Figs, [l] [5] 
|6) and [To] and the same shading is used in these figures. 



compare to Figs. [T] and [6] The major features in the 
distributions of b and q persist despite a bigger number 
of random parameters defining the Hamiltonians, more 
complex geometry of the two-level models, and a more 
chaotic resulting dynamics. There is a peak in the dis- 
tribution of the fractional collectivity b near 1 indicat- 
ing a sizable number of collective cases. The distribu- 
tion of the quadrupole moment for the collective real- 
izations (shaded in red) has a well-defined peak on the 
oblate side. The non-collective realizations appear to 
have quadrupole moment distribution centered at zero 
(shaded in blue). 

For systems with exact particle-hole symmetry the 
quadrupole moment for particles is equal in magnitude 
and opposite in sign to that of holes. Moreover, prop- 
erties such as excitation energies, spins of states, and 
transition rates, are exactly equal for particle-hole con- 
jugated systems. The particle-to-hole transformation for 
any two-body Hamiltonian amounts to the same Hamil- 
tonian for holes but with an additional one-body term. 
Thus, the symmetry is not exact in a two-level model 
space. Nevertheless in the TBRE, where two-body ma- 
trix elements are selected symmetrically about zero, the 
one-body term averages to zero. Therefore the results 
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0.2 0.4 0.6 0.8 1 -1 



b q 

Figure 12: (Color online) (13/2+, 13/2+)^^ (a) The distribu- 
tion of the fractional collectivity 6. (b) The distribution of the 
intrinsic quadrupole moment q. The system is particle-hole 
conjugated to that in Fig. |10[ The percentage of samples 
with the 0gs,2i sequence is 8.1% which includes 4.8% of col- 
lective and 1.7% of non-collective. This figure is analogous 
to Figs. [Tj [S] [6l [To] and |ll| and the same shading is used in 
these figures. 
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b q 

Figure 13: (Color online) (O/7/2, lp3/2)*- (a) The distribution 
of the fractional collectivity b. (b) The distribution of the 
intrinsic quadrupole moment q. The solid black line outlines 
the probability distribution for 31% of realizations with the 
Og.s. state followed by the 2i first excited state, both states 
with isospin T = 0. The 8.8% of realizations are collective and 
the 12.8% are non-collective. This figure is analogous to Figs. 
[1] [5] [6] [T0p2l and the same shading is used in these figures. 



in Figs. [12] and 10 for particle-hole conjugates systems 
(13/2+,13/2+)22^nd (13/2+,13/2+)6 are nearly sym- 
metric. The main difference is that the Hamiltonian 
for holes contains random single-particle energies which 
leads to a different fraction of collective realizations in 
the ensembles. 



V. REALISTIC VALENCE SPACE 

The schematic models discussed in the previous sec- 
tions all seem to possess the rotational low-lying spec- 
trum which is an evidence of the intrinsic deformation. 
However, to what extend they reflect the dynamics of 
realistic nuclei remains a question. The oblate intrinsic 
deformation observed in our models seems to be incon- 
sistent with the prolate dominance in real nuclei (see dis- 
cussion in Sec. VI I, moreover, the semi-magic nuclei with 



In Fig. 



one type of valence nucleons are generally not deformed. 
To attend to these issues we examine a realistic valence 
space consisting of the O/7/2 and single-particle lev- 
els, allowing for both protons and neutrons. The matrix 
element of the quadrupole operator for this model are 
constructed using the harmonic oscillator single-particle 
wave functions, we use the same effective charge for both 
types of nucleons. The multipole operator in this form 
facilitates comparison with the SU(3) group. 

we present our results for the (O/7/2, 
8 nucleons: 4 protons and 4 neutrons. This 



In Fig 
system wi 



corresponds to the configuration space of ''^Cr nucleus. 
In Fig. [l3ja), where the fractional collectivity h is shown, 
a noticeable peak that corresponds to collective realiza- 
tions is observed. The distribution of the quadrupole 
moment in Fig. 13 h) shows prolate and oblate peaks. 



The peaks are especially clear for the collective realiza- 
tions (shaded in uniform red) . The non-collective cases in 
Fig. 13 b) are distributed around g = (shaded in blue 
pattern). In agreement with the results in Ref. j7], hi 
this TBRE the prolate intrinsic shape is more probable, 
as evident from a bigger prolate peak. 



14 we focus on the 8.8% of realizations that 
The quadrupole moments in Fig. 13 b) 



are collective 

are further separated into prolate q > 0.7 and oblate 
q < —0.7 cases as shown in the inset of Fig. [T4j The 
same shading is used in the main figure showing the dis- 
tribution of the relative transition strength s. The max- 
imum possible value s = 1 is reached when the ground 
state wave function of the randomly selected Hamilto- 
nian coincides with that of the QQ Hamiltonian. From 
the summary in Tab. [Tjone finds that the J = 0, T = 
ground state of the QQ Hamiltonian, for which s — \, 
is prolate in this valence space. Indeed, the distribution 
of prolate realizations is peaked at around s = 0.8, while 
the oblate shapes have s near s = 0.6. 

The distributions of B^2 and i?42 for collective realiza- 



tions are shown in Fig. 15 This figure can be compared 
to Fig. |3l In both fi gures we use the same shading to 
separate the prolate and oblate collective cases. In con- 
trast to Fig. [3ja), both prolate and oblate realizations 



in Fig. 15 'a) have a band structure with the deexcita- 
tion ratio B^2 that is consistent with the rotational value 
10/7. This ensemble, based on the more realistic model 
space, appears to have an energy spectrum that is closer 
to the rotational spectrum. The distribution of i?42 in 
Fig. [T5|b) is broad, but it has a peak around the rotor 
value of 10/3. 

As concluded in Ref. [7J, realizations with rotational 
features appear in random ensembles due to correlated 
interaction matrix elements. Similarly to the single j 
level model, it is natural to attribute this collectivity 
to the QQ component in the Hamiltonian. The over- 
lap X between the ground state wave functions of the 
two-body random ensemble |Ogs(TBRE)) and the fixed 
ground state wave function of the QQ Hamiltonian is 
defined as 



= |(0<„(TBRE)|Og,(QQ))|^ 



(17) 



Fig. [16] shows the distribution of the overlap x in the 
(0/7/25 1^3/2)'^ model. A similar approach has been used 
in investigations of pairing coherence in random ensem- 




Figure 14: (Color online) (O/7/2, 1^3/2)**- The distribution of 
the relative transition strength s for the collective realizations 
(shaded with uniform red in Fig. 13 1. The 3.6% of prolate 



cases and 1.0% of oblate are identified with shades of color 
and pattern (see the inset). This figure is analogous to Fig. 
[2j and the same shading is used as in Figs. |2|4[ 
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Figure 15: (Color online) (O/7/2, lP3/2)*- (a) The distribution 
of the deexcitation ratio B42 defined in Eq. Q. (b) The 
distribution of the excitation energy ratio 7?42 defined in Eq. 
Q. Collective realization discussed in Fig. 



13 



are selected 

and, in addition, we require that the second excited state has 
spin 4. The fraction of such cases is 4.2%, with 2.4% being 
prolate and 0.6% being oblate, they are shaded separately 
with the same patterns as in Fig. |14[ The values for ,642 and 
7?42 from the QQ Hamiltonian listed in Tab |I] are shown with 
the vertical grid lines. This figure is analogous to Fig. [3] and 
the same shading is used as in Figs. |2|4|and|14| 



bles, see review in Ref. [2 . We select 56.3% of real- 
izations where the ground state quantum numbers are 
Jgs = and Tgs = 0, the ground state of the QQ 
Hamiltonian has the same spin and isospin. The dis- 
tribution of X shown in Fig. [16] is compared with the 
Porter-Thomas distribution. The latter emerges for 
uncorrelated wave functions in the 126-dimensional space 
spanned by the J = 0, T = wave functions. As shown in 
Fig. [16] the Porter-Thomas distribution drops abruptly, 
thus predicting that cases with large x are extremely un- 
likely. According to the Porter-Thomas distribution the 
probability to find a; > 0.1 is only 0.03% whereas in the 
TBRE X > 0.1 in 18.8% of random realizations. To em- 
phasize the relation between the collective structure and 



Figure 16: (Color online) (O/7/2, 1 ^3/ 2)^- The distribution 
of the overlap x defined in Eq. ( |17[ ). The results for all 
Jgs ~ 0, Tgs — states are unshaded; the fraction of such 
realizations is 56.3%. Collective realization that in addition 
to the Ogs state have the T = 0, 2i first excited state and b > 
0.7 are shaded (their fraction is 8.8%). Solid line shows the 
Porter-Thomas distribution, which is expected for the overlap 
between uncorrelated states. 



the large QQ component of the wave function we show 
in Fig. 16 the histogram for collective realization (with 

It is clear that the col- 



states Ogs and 2i and b > 0.7). 
lective transitions and rotational structure emerge when 
the component of the wave functions that corresponds to 
the eigenstate of the QQ Hamiltonian is large. 



VI. SUMMARY 

Our studies show that a collective behavior that re- 
sembles realistic is quite likely to be present in the en- 
semble with two-body random interactions. This behav- 
ior appears to emerge due to the quadrupole-quadrupole 
interaction component in the Hamiltonian. This com- 
ponent, as well as some higher multipoles can establish 
some noticeable coherence despite the overall many-body 
randomness and complexity. Similarly to the moment-of- 
inertia-like term (that is responsible for the ground 
state configurations with the maximum possible spin) 
the QQ component, while not a constant of motion, is 
dynamically prevailing. Let us list the supporting argu- 
ments: 

• The fraction of random realizations that are 
quadrupole-coUective is extremely large as com- 
pared to the statistically expected number. 

• In the two-body random ensemble, the quadrupole 
collectivity displayed by the transition rates disap- 
pears when the QQ component in the interaction 
is removed, see Fig. [9] 
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b q B42 R42 
(19/2)"^ 0.97 -0.979 1.42 3.31 

(19/2)* 0.95 -0.969 1.43 3.27 

(13/2+, 13/2+)^ 0.98 -0.977 1.41 3.29 
(13/2+,13/2-)® 0.98 -0.977 1.41 3.29 
(0/7/2,1^3/2)'* 0.97 0.996 1.35 3.27 

Table I: Characteristics of the QQ Hamiltonian. Listed in the 
table are the values of the fractional collectivity 6, quadrupole 
moment q, ratios of the transition rates and the ratios of 
excitation energies i?42. The models are the same as those 
considered in our study of the TBRE. 

• From investigations in Fig. [16] as well as indirectly 
form Figs. |2][7] and[T4] it follows that the collective 
states in the TBRE have structure similar to that 
of the QQ Hamiltonian eigenstates. 

• In order to examine the shape and other quantita- 
tive characteristics of the deformed mean-field we 
turn to the QQ Hamiltonian, for which the geome- 
try of the configuration space is the only parameter. 
The values of the quadrupole moments, transition 
rates, and level spacings for the models discussed in 
this paper are summarized in Tab. [l] In all cases the 
QQ Hamiltonian has a low-lying rotational spec- 
trum. The type of the quadrupole deformation and 
most of the other quantitative measures in Tab. [T] 
are consistent with those observed in the TBRE. 
This again suggests that the collective features seen 
in the TBRE arise from the coherent QQ compo- 
nent. 

Practically all deformed nuclei in nature are known 
to have a prolate ground state shape. This prolate domi- 
nance has been widely discussed in the literature |20fl25] . 
An effort to pinpoint the origin of the phenomenon us- 
ing the shell model approach with random interactions is 
presented in Ref. While in this work we do not ex- 
plicitly pursue the question of prolate dominance, we are 
compelled to comment on the issue from the standpoint 



of our findings. Our studies fully confirm the results in 
Ref. |7]. However, conclusions supporting the prolate 
dominance are difficult to draw, instead we offer several 
observations. 

First, the quadrupole collectivity seen in the TBRE 
is due to the QQ component in the Hamiltonian. This 
interaction and the geometry of the valence space deter- 
mine the deformation type. Thus, some questions of the 
shape systematics can be addressed by considering the 
QQ Hamiltonian and without invoking random interac- 
tions. 

Second, the shape is determined by the valence config- 
uration and by the positions of the single particle levels. 
The role of the single-particle level structure discussed by 
Hamamoto in Ref. |23| is possible to pinpoint using the 
TBRE as well as using analytic models, e.g. the seniority 
model and Elliot's SU(3) model [21]. 

Third, due to particle-hole symmetry, which does not 
need to be exact, the number of prolate and oblate con- 
figurations is approximately the same within a given va- 
lence space. The deviations from this symmetry affect 
only a few mid-shell systems where the two shapes com- 
pete. The effect of the particle-hole symmetry is seen in 
our results in Figs. [I0]and[l2] 

To conclude, in this work we examined the quadrupole 
collectivity that emerges in systems with two-body ran- 
dom interactions. A low-lying spectrum, characteristic 
of a rigid rotor, is commonly observed. The transi- 
tion B{E2,0gs — > 2i), the quadrupole moment of the 
first 2i state, and the deexcitation ratio B{E2,Ai — >■ 
2i)/ B{E2,2i — >■ Ogs) are all consistent with that of the 
deformed rotor. A weak triaxiality is also identified. The 
coherent dynamical role of the quadrupole-quadrupole in- 
teraction component is established as a source of this be- 
havior. 
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